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Abstract 

Q^ ■ We provide a new, independent, and analytic estimate of the lowest glueball mass, 

. and we found it at 1661 MeV within a relativistic quantum-field model based on ana- 

I lytic confinement. The conventional mesons and the weak decay constants are described 

5_i ' to extend the consideration. For the spectra of two-gluon and two-quark bound states 

^ ■ we solve the ladder Bethe-Salpeter equation. By using a minimal set of model param- 

eters (the quark masses, the coupling constant, and the confinement scale) we obtain 
' numerical results which are in reasonable agreement with experimental evidence in 

the wide range of energy scale. The model serves a reasonable framework to describe 
' r"] ', simultaneously different sectors in low-energy particle physics. 
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. 1 Introduction 

\D _ 

r^j- ■ Confinement and dynamical symmetry breaking are two crucial features of QCD, although 
Q . they correspond to different energy scales [1, 2]. Confinement is an explanation of the physics 
^ ! phenomenon that color charged particles are not observed; the quarks are confined with 
O I other quarks by the strong interaction to form bound states so that the net color is neutral. 
^ ' However, there is no analytic proof that QCD should be color confining and the reasons 
for quark confinement may be somewhat complicated. There exist different suggestions 
about the origin of confinement, some dating back to the early eighties (e.g., [3, 4]) and 
some more recent based on the Wilson loop techniques [5] , string theory quantized in higher 
5^ I dimensions [6], and lattice Monte Carlo simulations (e.g., [7]), etc. It may be supposed that 
the confinement is not obligatory connected with the strong-coupling regime, but it may be 
induced by the nontrivial background fields. One of the earliest suggestion in this direction 
is the analytic confinement (AC) based on the assumption that the QCD vacuum is realized 
by the self-dual vacuum gluon fields which are stable versus local quantum fluctuations and 
related to the confinement and chiral symmetry breaking [3] . This vacuum gluon field could 
serve as the true minimum of the QCD effective potential [8]. Particularly, it has been 
shown that the vacuum of the quark-gluon system has the minimum at the nonzero self-dual 
homogenous background field with constant strength and the quark and gluon propagators 
in the background gluon field represent entire analytic functions on the complex momentum 
plan [9]. However, direct use of these propagators for low-energy particle physics problems 
encounters complex formulae and cumbersome calculations. 

We are far from understanding how QCD works at longer distances. The well-established 
conventional perturbation theory cannot be used at low energy, where the most interesting 
and novel behavior is expected [10]. The calculations of hadron mass characteristics on the 
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level of experimental data precision still remain among the unsolved problems in QCD due to 
some technical and conceptual difficulties related with the color confinement and spontaneous 
chiral symmetry breaking. In such a case, it is useful to investigate the corresponding 
low-energy effective theories instead of tackling the fundamental theory itself. Although 
lattice gauge theories are the way to describe effects in the strong-coupling regime, other 
methods can be applied for some problems not yet feasible with lattice techniques. So 
data interpretations and calculations of hadron characteristics are frequently carried out 
with the help of phenomenological models. Different nonperturbative approaches have been 
proposed to deal with the long distance properties of QCD, such as chiral perturbation 
theory [11], QCD sum rule [12], heavy quark effective theory [13], etc. Along outstanding 
advantages these approaches have obvious shortcomings. Particularly, rigorous lattice QCD 
simulations [14] suffer from lattice artifacts and uncertainties and cannot yet give a reliable 
result in the low-energy hadronization region. The coupled Schwinger-Dyson equation (SDE) 
is a continuum method without IR and UV cutoffs and describes successfully the QCD 
vacuum and the long distance properties of strong interactions such as confinement and 
chiral symmetry breaking (e.g., [15]). However, an infinite series of equations requires to 
make truncations which are gauge dependent. The Bethe-Salpeter equation (BSE) is an 
important tool for studying the relativistic two-particle bound state problem in a field theory 
framework [16]. The BS amplitude in Minkowski space is singular and therefore, it is usually 
solved in Euclidean space to find the binding energy. The solution of the BSE allows to obtain 
useful information about the understructure of the hadrons and thus serves a powerful test 
for the quark theory of the mesons. Numerical calculations indicate that the ladder BSE 
with phenomenological potential models can give satisfactory results (for a review, see [17]). 

It represents a certain interest to combine the AC conception and the BSE method within 
a phenomenological model and to investigate some low-energy physics problems by using the 
path-integral approach. Particularly, it is shown that a "toy" model of interacting scalar 
"quarks" and "gluons" with AC could result in qualitatively reasonable description of the 
two- and three-particle bound states [18] and obtained analytic solutions to the ladder BSE 
lead to the Regge behaviors of meson spectra [19]. 

Below we consider a more realistic model introduced in [20] by taking into account the 
spin, color and fiavor degrees of constituents. This model was further modified in [21], 
applied to leptonic decay constants in [22], and used to simultaneously compute meson 
masses and estimate the mass of the lowest- lying glueball in [23, 24]. Here the aim is 
to collect all necessary formulae, explain the method in detail, and show that the correct 
symmetry structure of the quark-gluon interaction in the confinement region refiected in 
simple forms of the quark and gluon propagators can result in quantitatively reasonable 
estimates of physical characteristics in low-energy particle physics. In doing so, we build 
a model describing hadrons as relativistic bound states of quarks and gluons and calculate 
with reasonable accuracy the hadron important characteristics such as the lowest glueball 
mass, mass spectra of conventional mesons, and the decay constants of light mesons. 

The paper is organized as follows. In Sec. 11 we describe the main structure and specific 
features of the model. Analytic formulae for the meson spectra and weak decay constants are 
derived and numerical results on the vector and pseudoscalar meson masses and constants 
/tt and /k are evaluated in Sec. HI. The formation of a two-gluon bound state, the analytic 
expression for the lowest glueball mass and its numerical value are represented in Sec. IV. 
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2 The Model 



Because of the complexity of QCD, it is often prudent to examine simpler systems exhibiting 
similar characteristics first. Consider a simple relativistic quantum-field model of quark- 
gluon interaction assuming that the AC takes place. The model Lagrangian reads [23]: 

>C = {f^. - gf^^'^A^.A^f + E {q} [7„a" - + gTZ,A^^f ,}) , (1) 

where A^ - gluon adjoint representation (a = {!,... ,4}); = d^^A^ — d^A^; /^^*^ - 
the SUc{3) group structure constant {{A,B,C} = {1, ...,8}); - quark spinor of flavor / 
with color a — {1, 2, 3} and mass m/; g - the coupling strength, Pgr = ijat'^] and - the 
Gell-Mann matrices. 

Consider the partition function 

Z{g)^ jjVqVq JvAexpl^- JdxC[q,q,A]^ , Z(0) = 1 . (2) 

We allow that the couphng remains of order 1 (i.e., ag — 5'^/47r ~ 1) in the hadronization 
region. Then, the consideration may be restricted within the ladder approximation sufficient 
to estimate the spectra of two-quark and two-gluon bound states with reasonable accuracy 
[21, 23]. The path integrals defining the leading-order contributions to the two-quark and 
two-gluon bound states read: 

= JJVqDqexpi^-(qS-\) + ^ {(qrAq)(qrAq)) , (3) 

Zaa = (exp{-| > {(•))d = e-'^^^^-'^\.) . (4) 

The Green's functions in QCD are tightly connected to confinement and are ingredients 
for hadron phenomenology. The structure of the QCD vacuum is not well estabhshed and 
one may encounter difficulties by defining the explicit quark and gluon propagator at the 
confinement scale. Obviously, the conventional Dirac and Klein-Gordon forms of the prop- 
agators cannot adequately describe confined quarks and gluons in the hadronization region. 
Any widely accepted and rigorous analytic solutions to these propagators are still missing. 
Besides, the currents and vertices used to describe the connection of quarks (and gluons) 
within hadrons cannot be purely local. And, the matrix elements of hadron processes are 
integrated characteristics of the propagators and vertices. Therefore, taking into account 
the correct global symmetry properties and their breaking, also by introducing additional 
physical parameters, may be more important than the working out in detail (e.g., [25]). 

Because of the complexity of explicit Green functions derived in [9] , we examine simpler 
propagators exhibiting similar characteristics. Consider the following quark and gluon (in 
Feynman gauge) propagators: 

^ ^^^\ 2A^j ' 

^i''(p) = 5^''^exp(-pV4A% (5) 
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where p = p^Jn and u!{z) = (1 + z'^/i)'^. The sign "±" in the quark propagator corresponds 
to the self- and antiself-dual modes of the background gluon fields. These propagators are 
entire analytic functions in Euclidean space and may serve simple and reasonable approxi- 
mations to the explicit propagators obtained in [9]. Note, the interaction of the quark spin 
with the background gluon field generates a singular behavior S±{p) ~ l/""^/ in the massless 
limit m/ — > 0. This corresponds to the zero-mode solution (the lowest Landau level) of the 
massless Dirac equation in the presence of external gluon background field and generates a 
nontrivial quark condensate 

f fl^n r 1 fiA"^ r 777.^ 1 



(27r)^ 



indicating the broken chiral symmetry as mj ^ 0. A mass splitting appears between vector 
and pseudoscalar mesons (My > Mp) consisting of the same quark content. 

Our model has a minimal number of parameters, namely, the coupling constant as, the 
scale of confinement A and the quark masses {ttIu^;, tti^, mc, m;,}. Hereby, we do not make a 
distinction of the masses of lightest quarks, so m„ = rrifi = rriud- 

Below we describe the main steps in our approach on the example of the quark-antiquark 
bound state [24]. 

We allocate the one-gluon exchange between colored biquark currents 

^2 = |- 5^ JJ dxidx2 (qf,{xi)i-fi,t%,{xi)) D^^{xi,X2) {qh{x2)ilut^ qh{x2)) ■ (6) 

/1/2 

The color-singlet combination is isolated: 

9 3 

We perform a Fierz transformation 

(^7m)'5^'(^7.) =Y.Cj-OjOj, J = {5, P, V, A, T} , 
J 

where Cj = {1, 1, 1/2, -1/2, 0} and Oj = {/, ^75, ^7^, 757^, ^[7^, 7i.]/2}. 

For systems consisting of quarks with different masses it is important to pass to the 
relative co-ordinates (x, y) in the center-of-masses system: 

xi^x + $,iy, X2^x-^2y, Ci = , i = l,2. 

Then, we rewrite (6) 

2g_ 

Jhh 

where 



^2 = ^ E ^^// dxdyjj^,^,{x,y) D{y) jj;,;,(x, y), (7) 



Jjhf2ix,y) = (g/i(x + ^ly) Oj qf^{x - 62/)) • 
Introduce a system of orthonormalized functions {Uq{x)}: 

jdx Uq{x) Uq,{x) = , E Uq{^) Uqiy) = - y) ■ 

■' Q 
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Expand the biquark nonlocal current on the basis 

Div) J]m,{x, y) = ^Diy) j dz 5{z - y)^^ j]^^^^{x, z) 

= ^[dz y/m^pQiy) ■ 4d{z)Uq{z) J]f,f,{x, z) . 
Q 

We define a vertice function Vqjix^y^ 

Qfii^) VQj{x,y) (Ihi^) = \l D{y)UQ{y)qf^{x + iiy) Oj q}^{x - i2y) 

and a colorless biquark current localized at the center of masses: 

Jm{x)^ Jdy {qf,{x)VQj{x,y)qf,{x)) , J^^) ^ Jm{x) , Af^iQJfM- 
Then, (7) can be rewritten as follows: 

We represent the exponential by using a Gaussian path integral 

by introducing auxiliary meson fields Bj\f{x). Then, 

Z,g^mVqDqe^p{-{qS-'q)+g{Bj^Jj^)}^^ . 

Now we can take explicit path integration over quark variables and obtain 

Zgg ^Z^ (exp {IV In [1 + g{BMVM)S]})j, , 

where Tr = TrcTr^ J2±] Trg and Tr^ are traces taken on color and spinor indices, correspond- 
ingly, while J2± implies the sum over self-dual and antiself-dual modes. 

3 Mesons 

In particle accelerators, scientists see "jets" of many color-neutral particles in detectors 
instead of seeing the individual quarks. This process is commonly called hadronization and 
is one of the least understood processes in particle physics. 

We introduce a hadronization ansatz and will identify B^{x) fields with mesons carrying 
quantum numbers jV. We isolate all quadratic field configurations (~ Bjj-) in the "kinetic" 
term and rewrite the partition function for mesons [21]: 

Z^ [ YlVBj^ exp j-i J2 (Bj^ ['^^■^' + nA/-Ar'] B^f,) - Wres[BM] , (8) 
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where the interaction between mesons is described by the residual part W^resl-BA/"] ^ ^{^m)- 
The leading-order term of the polarization operator is 

^ATM'izi - Z2) = jjdxdy Uj^{x) asX{zi - Z2, x, y) Uj^,{y) , (9) 

where the Fourier transform of the kernel reads 

asXjAp,x,y) = asjdze^^^XjAz^x^y) = ^ D{x)D{y) e^'^^^-y^ 

•Tr [Oj~S (k + 6p) Oj,S [k - . (10) 
We diagonalize the polarization kernel on the orthonormal basis {Uf^}: 

JdyXjj'ip, X, y) U^,{y) = X^{-p^) U^,{x) 

or, 

j jdxdyUM{,x)Xjj>{p,x,y)Uj^i{y) =5^^' Xj^^{-p^) 

that is equivalent to the solution of the corresponding ladder BSE. 

In relativistic quantum-field theory a stable bound state of n massive particles shows up 
as a pole in the S-matrix with a center of mass energy. Accordingly, the meson mass may 
be derived from the equation: 

1 + asX:^{Ml) = , -p^ = . (11) 

The following renormalization takes place: 

{Um[1 + asX^{-p')]UM) = {Um[1 + ocsX^iMl) + asX^{M%)\p' + Ml^]U^f) (12) 

^{UR\p^ + Mlf\UR), Xm{z)^ 



dz 

where the renormalized state function reads 



Ur{x) = V«XK^ • Um{x) . (13) 

The use of the path-integral technique leads to the following practical advantages over 
simply solving a BSE with one-boson exchange: 

(i) the vacuum functional may be written in alternative representations, either through 
original variables of quarks and gluons or, in terms of bound states, i.e., we obtain so-called 
"quark-hadron duality" , 

(ii) the BS kernel (10) is natively obtained in a symmetric form, 

(iii) the normalization of the operators of bound states is performed in the most simple 
way by keeping the condition X{Mj) > evident, 

(iv) after renormalization (12) the partition function of the system of Bj\f fields takes the 
conventional form with a kinetic term and interaction parts. 
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3.1 Pseudoscalar and vector meson ground states 

In the quark model {qf^qf.^) bound states are classified in J'^'-^ multiplets. For a pair with 
spin s = {0, 1} and angular momentum £ the parity is P = (—1)^+* and the total spin is 
|£ — s| < J < \^ + s\. Below we consider the meson ground states {£ — 0, = 0), the 
pseudoscalar (P : J^^ — 0~+) and vector (V : J^^ — 1 ) mesons, the most established 
sectors of hadron spectroscopy. 

We should derive the meson masses from Eq. (11). The polarization kernel Xui^P^) is 
real and symmetric that allows us to find a simple variational solution to this problem. For 
the ground state we choose a trial function [21, 24]: 



[/(a;, a) ~ y D{x) ■ exp 



aA X 



2^2 



jdx\U{x,a)f — 1 , a>0. 



(14) 



Substituting (14) into (11) the variational equation defining the masses of P and V 
mesons as follows: 



1 = -a^ • Aj(A, Mj,mi,m2) 

= exp < — — — > max 

37rmim2 I 2A^ I i/4<a<i/2 



(6a-l)(l-2a) 



• exp 



2A2 



4apj + ^ + a(2 - apj){^, - 6)') 



A' 



mim2 / \ / \n 

H [1 + Xj (mi) a; (ms)] 



A2 



(15) 



where Cj = {1, 1/2}, pj = {1, 1/2} and Xj = {1, -1} for J = {P, V]. 

Localization of the meson field at the center of masses of two quarks results in the 
following asymptotic properties. For mesons consisting of two very heavy quarks (mi = 
m2 = m ^ 1) we solve (15) and obtain the correct asymptotic behavior 



ej = 41n 



37r 



^32(7-4^)Cja, 



Note, the next-to- leading value Sj does not depend on any masses. Moreover, Sy > ep 
because the corresponding Fierz coefficients obey Cp = 1 > Cy = 1/2. The mass splitting 
My > Mp remains for "heavy-heavy" quarkonia. 

For a "heavy- light" quarkonium (mi ^ 1 , m2 ~ 1) we estimate the mass 



ml-ej, ej^ej{Mj). 



3.2 Weak decay constants 

An important quantity in the meson physics is the weak decay constant. The precise knowl- 
edge of its value provides great improvement in our understanding of various processes con- 
volving meson decays. For the pseudoscalar mesons the weak decay constant fp is defined 
by the following current-meson duality: 



^/pPM = (0|^A(0)|^i^(p)), 
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where J a is the axial vector part of the weak current and Ur{p) is the normahzed vector of 
state. 

We estimate 



V2g f dk 



32Aa,^2A(M2) (1 - 2ap)(6ap - 1) 



3 y(27r) 
32 

• exp 



37r^/^(mi + 1712 



ml — m\ 



1 + 2apf 
ap 



1 + 



1 + 2ap 



dp {ttii — 1712] 

1 + 2ap 17111712 



(16) 



where ap is the value of parameter a calculated for the given meson with mass Mp. 

Particularly, for an "asymmetric" meson containing an infinitely heavy quark (mi ^ 
777,2 ~ 1) we obtain the correct asymptotic behavior 

fp ~ l/^JrrTi 

due to the localization of the meson field at the center of two quark masses. 



3.3 Numerical results 

To calculate the meson masses we need to fix the model parameters. We determine the 
quark mass rUud and the coupling constant ag from equations: 

1 + Q;sAp(A, 138 MeV, m„<i, rriud) = , 1 + oisAy (A, 770, m„d, rriud) = (17) 

by fitting the well-established mesons 7r(138) and p(770) at different values of A. The remain- 
ing constituent quark masses ms,inc, and mj, are determined by fitting the known mesons 
X(495), J/* (3097), and T(9460) as follows: 

1 + asXp{A, 495, niud, rUs) = , 
1 -|- Q;sAy(A, 3097, mc, rUc) — , 
1 q;sAv(A, 9460, m^, m^) = . 

The dependencies of the estimated constituent quark masses on A are plotted in Fig. 1. 

The sharp drop of all quark mass curves in Fig.l may be shortly explained as follows. 
Note, two equations in Eqs. (17) mostly differ by meson masses in exponentials along 
different numerical factors Cj, pj and x.J- They have general solutions {mud , (^s} not for 
any A. Suppose, at fixed A = Aq they arc solvable. Then, for finite coupling ct^ the solution 
mud is obviously finite to obey both equations. However, for vanishing ag ^ the equations 
take the form ^ 

1 ~ ^ o • const{Ao, Mj, Pj) 

Kd 

and the solution for quark mass behaves m^d ~ ^/«s ~^ 0. This picture is observed in Fig.l. 

By using these quark masses and coupling constant we can estimate other meson masses 
in dependence on A and some results are shown in Fig. 2. 

To fix the value of parameter A we calculate the weak decay constants and fx to 
compare with experimental data. Note, these constants considerably depend on A (see Fig. 
3) that allow us to fix it unambiguously at A = 416.4 MeV. 
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Figure 1: Solutions for constituent quark masses versa the confinement scale value A. 

The final set of model parameters are fixed as follows: 

as = 1.5023 , A = 416.4 MeV , m^^ = 206.9 MeV , 
rus = 323.6 MeV , = 1453.8 MeV , nib = 4698.9 MeV . (18) 

With these parameters we have estimated the pseudoscalar and vector meson masses shown 
in Table 1 and compared with experimental data [26]. The relative error of our estimate 
does not exceed 3.5 percent in the whole range of mass (from 0.14 GeV up to 9.5 GeV). 
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3097 
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547 


5,(5370) 


5551 
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5^(6286) 


6522 
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1022 


T(9460) 


9460 


Dsimo) 


1970 


r76(9300) 


9434 


D*(2010) 


1942 







Table 1: Estimated spectrum of conventional mesons (in units of MeV). 



There are mainly two schemes describing lj — $ and 77 — 77' mixings [26]. The octet-singlet 
scheme uses the mixing angle 9 between states {uu + dd—2ss)/^/6 and (uu+ dd+ ss) / ^/S. We 
use the quark-flavor based mixing scheme between states {uu + dd) / \pl and ss with mixing 
angle y?. These two schemes are equivalent to each other by6' = (/9 — 7r/2-|- arctan(l/\/2) 
when the S'f/(3) symmetry is perfect. Particularly, for "ideal" vector mixing the angle is 

= 90° or eif = 35.3°. 

With fixed parameters (18) we calculate a relatively heavy mass My(ss) = 1064 MeV of 
vector ss state. To obtain correct masses of to'(782) and $(1019) one needs a considerable 
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Figure 2: Solutions for some meson masses in dependence on the confinement scale value A. 

mixing to the light quark-antiquark state with mixing angle (fy ~ 73.2° which differs signif- 
icantly from the "ideal" value. By using the same parameters (18) we obtain a pseudoscalar 
ss state with mass Mp{ss) = 705 MeV. We cannot describe the physical mass of ri'{958) 
by any mixing to the light-quark pair and can only fit the correct mass Mp{ri) = 547 MeV 
at angle (fp — 58.5°. Our model fails to describe simultaneously the t] — t]' mixing. This 
problem obviously deserves a separate consideration. 

Note, the infrared behavior of effective (mass-dependent) QCD coupling a.^ is not well 
defined and needs to be more specified [27, 28, 29]. In the region below the r-lepton mass 
{Mr = 1.777 GeV) the strong-coupling value is expected between ^^(Mt-) ^ 0.34 [26] and 
the infrared fix point 0^(0) = 2.972 [30]. Our parameter = 1.5023 does not contradict 
this expectation because it is estimated to fit the vr meson mass, and so the corresponding 
energy scale is ~ 140 MeV. We keep this value for further calculations. 

The weak decay constants of light mesons are well established data and many groups 
(MILC [31], NPLQCD [32], HPQCD [33], etc.) have these with accuracy at the 2 percent 
level. Therefore, these values are often used to test any model in QCD. By substituting 
optimal values of {m^^rf, m^, a^. A} (18) into (16) we calculate 

U = 128.8 MeV , fx = 157.7 MeV . 

Our estimates are in agreement with the experimental data [34, 26]: 

/™ = 130.4 ± 0.04 ± 0.2 MeV , f^?^ = 155.5 ± 0.2 ± 0.8 ± 0.2 MeV . (19) 

Our model represents a reasonable framework to describe the conventional mesons, and 
the parameters are fixed. Below we can consider two-gluon bound states. 



4 Glueball Lowest State 

Because of the confinement, gluons are not observed, they may only come in bound states 
called glueballs. Glueballs are the most unusual particles predicted by the QCD but not 
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Figure 3: Weak decay constants depending on the confinement scale value A. 



found experimentally yet [35]. There are predictions expecting non-gg scalar objects, like 
glueballs and multiquark states in the mass range ~ 1500 — 1800 MeV [36, 37]. Experimen- 
tally the closest scalar resonances to this energy range are the /o(1500) and /o(1710) [38]. 
Some references favor the /o(1500) as the lightest scalar glueball [39], while others do so for 
the /o(1710) [40, 41]. Recent scalar hadron /o(1810) reported by the BES collaboration may 
also be a glueball candidate [42]. 

The study of glueballs currently deserves much interest from a theoretical point of view, 
either within the framework of effective models or lattice QCD. The glueball spectrum has 
been studied by using effective approaches like the QCD sum rules [43], Coulomb gauge QCD 
[44], and potential models (e.g., [45, 46]), etc. The potential models consider glueballs as 
bound states of two or more constituent gluons interacting via a phenomenological potential 
[45, 47, 48]. It should be noted that potential models have difficulties in reproducing all 
known lattice QCD data. Different string models are used for describing glueballs [49, 50], 
including combinations of string and potential approaches [46]. It has been shown that a 
proper inclusion of the helicity degrees of freedom can improve the compatibility between 
lattice QCD and potential models [51]. 

An important theoretical achievement in this field has been the prediction and computa- 
tion of the glueball spectrum in lattice QCD simulations [52, 53] . Recent lattice calculations, 
QCD sum rules, "tube" and constituent glue models predict that the lightest glueball has 
the quantum numbers of scalar {J^'-' = 0++) and tensor (2"'"+) states [54]. Gluodynamics 
has been extensively investigated within quenched lattice QCD simulations and the lightest 
glueball is found a scalar object with a mass of ~ 1.66 ± 0.05 GeV [55]. A use of much finer 
isotropic lattices resulted in a value 1.475 GeV [53]. Recently, an improved quenched lattice 
calculation of the glueball spectrum at the infinite volume and continuum limits based on 
much larger and finer lattices have been carried out and the scalar glueball mass is calculated 
to be 1710 ± 50 ± 80 MeV [56]. 

Two-gluon bound states are the most studied purely gluonic systems in the literature, 
because when the spin-orbital interaction is ignored = 0), only scalar and tensor states are 
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allowed. Particularly, the lightest glueballs with positive charge parity can be successfully 
modeled by a two-gluon system in which the constituent gluons are massless helicity-one 
particles [57]. 

Below we consider a two-gluon scalar bound state. We isolate the color-singlet term in 
the bi-gluon current in (4) by using the known relations 

jABE jA'B'E _ 2 ^^AA' ^BB' _ ^AB' ^BA'^ _|_ ^AA'E^BB'E _ ^AB'E^BA'E 

The second-order matrix element containing a color-singlet two-gluon current reads [23] 

Laa= ^JJ dxdy {ji^ix,y)jZf(x,y) - jPix,y)J^^^ (x^y)) 

■ S'^^'W^^^ix, y) - Si'^'W^^^x, y) - S^^^'W^^x, y) + y) 
where 

J%^{x,y)^A^^{x)A<'Ay), 

W,.{x, y) - ^^^(^ -y) = 6>'^W{x-y) + ... , W{z) = ^ e"^' . 

This part consists of spin-zero (scalar) and spin-two (tensor) components. Below we consider 
the scalar component: 

^AA^ ^ jj dXidX2 J {xi,X2)W{xi- X2) J {xi,X2) , J{xi,X2) = JI^i^{xi,X2) . 

By introducing the relative coordinates {xi ^ x + y /2 , X2^ x — y/2) we rewrite 

2 

Laa^ |- jjdxdy J{x, y)W{y)J{x, y) . (20) 

One can see that the matrix element (20) is similar to (7) by the very construction. By 
omitting details of intermediate calculations (similar to those represented in the previous 
section) we rewrite the partition function in terms of auxiliary field B{x) as follows: 

Zaa^Zg = j VBe^v[-\ {bG-'b) + Li[B] 
where Li[B] ~ 0{B^) and the BS kernel is 

G-\x-y)^5{x-y)-^-Il{x-y), 

Tl{z) = JJdtds ^n(t) D + z^D - .Jw{s) Ur^is) . 

The hadronization ansatz allows us to identify B with scalar glueball field. To find the 
glucball mass we should diagonalize the Bethe-Salpeter kernel ll{z). The glueball mass Mq 
is defined from equation [24]: 



1-^ 



J dz e''P U{z) = , p^ = -Ml . (21) 
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For the lightest ground-state scalar glueball choose a Gaussian wave function: 



U{x) = - e-"^" 



dx\U{x)\' 



c> 0. 



Then, we derive (21) as follows: 

1 



a. 



exp 



4A2 



Oicrit 



37r(3 + 2V2)^ 



The final analytic result for the lowest-state glueball mass reads 

Mg = 2A 



a. 



(22) 



The solution Mq > exists for any < Ocrit ~ 80.041. 

Note, the scalar glueball mass depends linearly on the confinement scale A and the 
scaled mass Mq/A depends only on coupling as (see Fig. 4). Particularly, if we take values 
A ~ Aqcd ~ 360 MeV and a, ~ a,(M^) = 0.343, then we estimate Mg ^ 1710 MeV. 

However, our purpose is to describe simultaneously different sectors of low-energy particle 
physics. Accordingly, with values = 1.5023 and A = 416.4 MeV determined by fitting the 
meson masses and weak decay constants, we calculate the scalar glueball mass as follows: 



Mg = 1661 MeV . 



(23) 



Our estimate (23) is in reasonable agreement with other predictions expecting the lightest 
glueball located in the scalar channel in the mass range ~ 1500 — 1800 MeV [36, 43, 53, 58]. 
The often referred quenched QCD calculations predict 1750 ± 50 ± 80 MeV for the mass of 
the lightest glueball [52] . The recent quenched lattice estimate with improved lattice spacing 
favors a scalar glueball mass Mg = 1710 ± 50 ± 58 MeV [56]. 




Figure 4: Evolution of the lowest-state glueball mass scaled to A with the coupling as 
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Another important property of the scalar glueball is its size, the "radius" which should 
depend somehow on the glueball mass. We estimate the glueball size by using the "efTective 
potential" W{y) (20) connecting two scalar gluon currents. The glueball radius may be 
roughly estimated as follows 



Jd^x x"^ W{x) _ V2 1 
\ Jd^x W{x) - X ^ 295 MeV ""^^ " 



This means that the dominant forces responsible for binding gluons must be provided by 
medium-sized vacuum fluctuations of correlation length ~ 0.7 fm. Consequently, typical 
energy-momentum transfers inside a scalar glueball occur at the QCD scale ~ 360 MeV, 
rather than at the chiral symmetry breaking scale ~ 1 GeV (or, ~ 5fm). 
Prom (22) and (24) we deduce that 



In 



/ O^crit \ 



1/2 



5.64. 



This value may be compared with the prediction (r^ ■ Mq = 4.16 ± 0.15) of quenched QCD 
calculations [52, 56]. A study of the glueball properties at finite temperature using SU(3) 
lattice QCD at the quenched level with the anisotropic lattice imposes restrictions on the 
glueball parameters at zero temperature: 0.37 fm < < 0.57 fm and Mq ^ 1.49 GeV 
[59] . The nonprincipal differences of quenched lattice QCD data from our estimates may be 
explained by the presence of quarks (our parameters have been fixed by fitting two-quark 
bound states) in our model. 

A method of analysis of correlation functions in QCD is to calculate the corresponding 
condensates. The value of the correlation function dictates the values of the condensates. 
We calculate the lowest nonvanishing gluon condensate in the leading-order (ladder) approx- 
imation: 

y2Tr(F^^F7) = 8Nc7rasA.^ jd'^zW{z) = 67ra,A^ !v 0.8 GeV^ 
which is the same order of magnitude with the reference value [60] 

g^Tr {Gi.M^"') ^ 0.5 GeV^ . 



In conclusion, the suggested model in its simple form is far from real QCD. However, our 
aim is to demonstrate that global properties of the lowest glueball state and conventional 
mesons may be explained in a simple way in the framework of a simple relativistic quantum- 
field model of quark- gluon interaction based on analytic confinement. Our guess about the 
symmetry structure of the quark-gluon interaction in the confinement region has been tested 
and the use of simple forms of propagators has resulted in quantitatively reasonable estimates 
in different sectors of the low-energy particle physics. The consideration can be extended to 
other problems in hadron physics. 

The author thanks G.V. Efimov, E. Klempt and V. Mathiew for valuable remarks and 
useful suggestions. 
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